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We present an experimental system to study the Bose polaron by immersion of single, well- 
controllable neutral Cs impurities into a Rb Bose-Einstein condensate (BEG). We show that, by 
proper optical traps, independent control over impurity and BEG allows for precision relative posi¬ 
tioning of the two sub-systems as well as for independent read-out. We furthermore estimate that 
measuring the polaron binding energy of Frohlich-type Bose polarons in the low and intermediate 
coupling regime is feasible with our experimental constraints and limitations discussed. 


I. INTRODUCTION 

The immersion of single, controllable atoms into a Bose-Einstein condensate (BEG) realizes a paradigm of quantum 
physics - individual quantum objects interacting coherently with a single or few mode bath. This system allows to 
experimentally address various questions of quantum engineering, including local, non-demolition measurement of 
a quantum many-body system [T]; cooling of qubits while preserving internal state coherence Hi; or engineering 
bath-mediated, long-range interaction between two or more impurities [3]. An impurity strongly interacting with 
the quantum gas will loose its single particle properties and it is rather described in terms of quasi particles, which 
are known as Bose polarons m- Particularly for strong interaction, such systems have been predicted to show 
remarkable properties such as self-trapping or polaron clustering [S]. Experimentally, impurities in BECs have 
been introduced, for example, as many atoms of different internal state mm or different atomic species [IT1[T1], or 
as individual ions [131E! or electrons m- 

The system we report on here considers a single, neutral impurity in an ultracold gas, where the effect of impurity- 
impurity interactions, either direct or mediated by the bath, can be neglected. The corresponding Hamiltonian thus 
reflects an extremely imbalanced mixture and can be written as ISdS] 
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The first row of Eq. (fy represents the BEG part of the Hamiltonian, where a|^ (a^.) creates (annihilates) a boson of 
mass to_b, momentuniK, and energy dispersion ek = / {2mB)- The interaction of bosons within the BEG is given 

by the contact s-wave interaction with Fourier transform yBB(q). The impurity atom of momentum pi, mass m/ and 
density pi as well as its interaction with the BEG via the potential Vib are described by the second row in Eq. ([^. 


A. The Frohlich polaron 

Originally, the polaron concept was developed for condensed matter systems to describe electrons moving in a crystal 
lattice. The interaction between a moving electron and lattice of the ion cores forms a propagating quasi-particle 
called polaron, comprising electron and surrounding phonon cloud. Specifically the effective mass of the polaron as 
well as its energy could strongly differ from the bare electron’s values depending on the interaction with the crystal 
[13 HH]. In the limit of very small electronic wave vector, the underlying crystal structure can be neglected, and the 
crystal can be described as a continuously polarizable medium. In this case, the model to describe polarons is the 
well-known Erohlich Hamiltonian [19] 
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given by the sum of the kinetic energy of the impurity particle with mass m/, the energy of the phonons in the 
medium with dispersion Wk, and the interaction energy with coupling constant Vk arising between the two of them. 
The coupling strength is essentially determined by the dielectric constant of the crystal and can be quantified by a 
single, dimensionless coupling constant a. The regime of small a is characterized by weak and intermediate polaron 
coupling while for values of a 3> 1 the strong coupling regime is realized. Simple perturbation theory predicts the 
crossover around a value of a « 6. Experimentally, Frohlich polarons have been observed in condensed matter systems 
for small and moderate interaction strengths a, however, the regime of so-called strong coupled polarons with a ^ 1 
has so far not been experimentally investigated mMi- 

Recently, the quantum gas Bose-polaron was subject of intense theoretical work. It was shown that the description 
of a Bose-polaron Eq. Q can be directly mapped onto the Frohlich Hamiltonian Eq. ([2 1 via a Bogoliubov transform [5] 
and that this model could be accessed experimentally dill]. This analogy holds as long as the interaction between 
two BEC excitations can be neglected. The resulting Hamiltonian has the same operator structure as Eq. Q, 
where the crystal (optical) phonons are replaced by elementary Bogoliubov excitations of the BEC with dispersion 
Eq = hvag^yi + 9 )^/ 2 , where Vg = h/(\/2m^) is the sound velocity in the BEC, ^ = l/V87rnBEC cibb is the healing 

length, riBEC is the BEC atomic density, and obb is the s-wave scattering length between two atoms in the BEC. 
The polaronic coupling parameter in the dilute quantum gas scenario PE] 
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is given by the s-wave scattering lengths between a Boson of the quantum gas with the impurity (oib) or with other 
quantum gas atoms (obb), and the condensate healing length For the quantum gas Bose-polaron, these parameters 
can be varied over a large range of values in the weakly interacting regime and the strong coupling regime by means 
of Feshbach resonances. Moreover, the characteristic properties of the emerging polaron which are the binding energy 
Ep and its effective mass rup can be inferred by, e.g., radio frequency spectroscopy and trap frequency measurements, 
respectively m- Importantly, also the Bose polaron beyond the Frohlich model can be experimentally realized 
yielding access to a rich and highly controllable model system of impurity physics in quantum fluids. 


II. EXPERIMENTAL REALIZATION 

Our experimental approach to realizing the Bose polaron aims at immersing single or few neutral Caesium (^^^Cs) 
atoms into a Rubidium (®’^Rb) BEC. This combination of species features several advantages, facilitating the realiza¬ 
tion, control and characterization of the Bose polaron. First, due to the relatively high nuclear charge of Rb and Cs 
the fine structure splitting of both species is also relatively large and allows to tune dipole traps in between the two 
fine structure lines of the first excited P-level with moderate unwanted photon scattering [521 !23] • As a consequence 
the atoms of this element do not experience a dipole potential. In order to improve the control over both species inde¬ 
pendently, we employ this fact constructing a species-selective lattice allowing for trapping and controlled transport 
of impurity atoms, only. 

Second, for a dipole trap wavelength of A = 1064 nm, the trapping frequency w and thus the gravitational sag g/co'^ 
with g the gravitational acceleration, is equal to the percent level for the two species. 

Third, using few or even single impurity atoms has several advantages. For Cs representing the minority component, 
three-body losses limiting the lifetime of the polaron are due to Rb-Rb-Cs collisions rather than Cs-Cs-Rb collision, 
where the loss coefficient L 3 of the former is an order of magnitude smaller than for the latter [24]. While three-body 
losses will still be a limitation of the atomic lifetime, the lifetime of the polaron is expected to be significantly larger 
than the lifetime in balanced Rb-Cs mixtures. However, if the elastic interaction of multiple impurities is desired 
this can be controlled by adjusting the number of trapped impurities. Furthermore, fluorescence imaging of such a 
small number of atoms in an optical lattice allows for single-site resolved detection of the impurities with standard 
optical systems |25| . Moreover, only the dynamics of a single impurity allows to dynamically track the trajectory of 
a polaron, because crossing of the trajectories of multiple indistinguishable impurities would lead to an ambiguity. 

For the design of the experimental apparatus, additional considerations have to be made. In a combined system 
of quantum gas and single atoms, the respective ways to experimentally extract information from averages differ: 
For a quantum gas, a single realization yields an ensemble average of typically 10^ ... 10® atoms. For single atoms, 
in contrast, averages have to be formed as time averages of typically 10 ^ ... 10 ® repetitions for identical parameter 
values. For a combination, the statistics is clearly limited by probing the single impurity, while usually the time scale 
of a single experimental run is limited by the production of a BEC. 

In the following we first discuss the experimentally relevant parameter ranges and constraints for our system of 
single Cs impurities in a Rb BEC, before we turn to the presentation of our experimental apparatus. 
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FIG. 1. Bose polaron properties at the Feshbach resonance. 

a) shows the theoretical values of the interspecies scattering length aiB and the three-body loss coefficient K 3 for Rb-Rb-Cs 
scattering at the 181 G Feshbach resonance (dashed red). For our typical BEG density of 10^^ cm“^ the estimated coupling 
strength a and the corresponding polaron energy Ep is shown in b). The vertical blue shaded area indicates fulfillment of the 
condition j3 = r'p/yp > 5 (> 1 light shaded), c) The parameter e is derived in [7] and gives an estimate of how well the Bose 
polaron fits the Frohlich type. The horizontal green shaded area indicates e < 0.3 which means good correspondence to the 
Frohlich type polaron. 


The Rb-Cs Bose-Polaron 

For the experimental characterization of the Bose-Polaron, we focus on the binding energy = hvp in the following. 
This energy can be measured by radio or microwave spectroscopy, driving a Zeeman or hyperfine transition between 
two impurity states, where one state is interacting with the bosonic bath forming the polaron, whereas the other state 
is non-interacting [261E7]. The decay rate of the polaron 7 p implies a lower bound for the linewidth of the polaron 
spectroscopy peak and is dominated by three-body losses, i.e. molecule formation, occurring with rate 7 p = 
with the loss coefficient L 3 and the BEG density ubec- In order to clearly resolve the polaron peak at a frequency 
shift of Vp, the ratio 13 = Vpj^p should be significantly larger than 1 yielding a figure of merit for the determination 
of optimum experimental parameters. 

In Figure [n] we explore the regimes of Bose polarons realized with our typical experimental parameters discussed 
below, employing theoretical data of elastic scattering length oib and three-body loss coefficient L 3 around the 
interspecies Feshbach resonance at 181 G [55]. The range of polaronic coupling strengths a accessible with our 
experimental setup is determined by several parameters: the boson-boson scattering length obb ~ lOOoo |29j . with 
the Bohr radius oq = 53 pm, is given by the background scattering length of Rb and does not change within the range 
of magnetic field values considered here; BEG peak densities are on the order of riBEC.max = 10^"^ cm“^ which serves 
as a worst case approximation for riBEc; and the interspecies scattering length oib ~ 650 cq at zero magnetic field is 
tunable in the vicinity of the Feshbach resonance (301 ED. 

From this we calculate the polaron coupling constant a and polaronic binding energy Ep using a simple analytic 
expression for an impurity with infinite mass from |16j . see Figure [TT|)b). Furthermore, the polaronic decay rate 7 p is 
calculated from the BEG peak density and theoretical values for L 3 in vicinity of the Feshbach resonance [28]. For 
/3 > 1, we can reliably resolve polarons spectroscopically within the discussed parameter ranges, indicated by the 
vertical shaded regions in Figure [ll] 
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Furthermore, in Figure [TTj[c) we use the parameter e from [7] to indicate to which extend the Bose polaron realized 
can be described by means of the Frohlich model discussed above, where e ^ 1 corresponds to a good description. 
One can see that e lies well below 1 for the experimentally directly accessible range (blue). For reference, the value of 
e < 0.3 which is given in [7] is marked as horizontally shaded area in Figure [Tl| 

For observation of a strongly coupled Frohlich polaron, clearly, the parameter range has to be adjusted by, for 
example, employing Feshbach resonances with more favorable three-body loss properties or optical tuning methods 
[32] . by increasing the sensitivity of the spectroscopic detection for a given life-time, and by adapting the BEC 
density via changes of the trap geometry. However, a Bose-polaron simulating the Frohlich model in the weak and 
intermediate coupling regime is well accessible by immersing single neutral Cs atoms into a BEC with density of the 
order of 10^"^ cm“^. 


Rb Bose-Einstein Condensate 


In order to optimize the statistics of single atom probing, the experimental setup aims at a short BEC production 
time. This is realized by, first, a short initial laser cooling stage, where a 3D magneto-optical trap (MOT) is loaded 
from a 2D MOT in « 2.5 s and, second, evaporation in a steep optical dipole trap which is formed by a horizontal 
and a vertical beam within « 4 s. In order to avoid perturbation from cooling and trapping of single Cs atoms in a 
MOT, the Rb cloud is prepared in the magnetic field insensitive F = l,mF = 0 state while evaporating. We typically 
prepare a BEC with 2.5 x 10"^ atoms at a peak density of 1.2 x 10^^ cm“^ and a critical temperature of approximately 
120 nK. The BEC’s decay rate of 0.4 Hz is dominated by two- and three-body collisions at this density but decreases to 
0.08 Hz for a lower number of atoms. For technical details of the BEC production and state preparation, see Section 

m 


Single Atoms 


Single or few Cs atoms are captured in a high magnetic field gradient (« 250G/cm) MOT [S^, spatially overlapped 
with the Rb MOT, but operated at a different time in the preparation sequence. In contrast to the Rb MOT, the Cs 
MOT is loading atoms from the background, and the laser cooling beams have a smaller diameter and beam intensity, 
leading to an overall reduced loading rate. The fluorescence signal of such a MOT features discrete values, which can 
be assigned to a specific number of atoms in the MOT. A Poissonian distribution describes the loading statistics in 
the few atom regime, observed in a histogram of measured fluorescence signal, corresponding to the distribution of 
atom numbers in Figure |h] (b). The width of each histogram peak is ideally given by the shot noise of the atoms’ 
fluorescence light; practically, however, fluctuations of the MOT laser beam intensities and technical contributions 
such as readout noise additionally broaden the peaks. This limits the maximum countable atom number to roughly 
eight for our setup. Refer to Section m for technical details and a description of the fluorescence imaging system. 


In order to provide position control over the single Cs atoms independently from the Rb BEC trap, we apply a 
species selective optical conveyor belt lattice, formed by two counter propagating, linearly polarized laser beams with 
wavelength Aiat = 790.0180 nm and a waist of 29 /im. For this wavelength between the Rb D-Lines, the resulting 
potential cancels out for Rb atoms [H1I15], but at the same time the frequency is blue detuned for Cs, providing tight 
confinement along the lattice axis in the nodes of the standing wave with depths up to 7300 = 850 fiK x ke, with 

/2mcs the single photon recoil energy, k = 27r/Aiat, ke the Boltzmann constant, and mcs the mass of a 
Cs atom. While the lattice provides tight axial confinement for Cs atoms, it does not confine the atoms radially. We 
therefore superpose the lattice axis with one beam of the dipole trap (see Section HI) and obtain a maximum trap 
depth of 1.45 mK radially, resulting in trap frequencies of 27r x 4 kHz radially and 27r x 460 kHz axially. The lifetime 
of atoms in the lattice at full depth is limited to r « 1.24 s by phase fluctuations. If in addition an optical molasses 
is used to cool the atoms, the lifetime can be extended up to r = 71 s, limited by background pressure. 


An important quantity characterizing the lattice is the selectivity s = E^^/E^^ for a given intensity and wavelength. 
By performing Raman-Nath scattering on the BEC for various lattice wavelengths, we have identified the 

optimal wavelength Aiat.vac = (790.0180 ± 0.0017) nm and find a selectivity of 5000 : 1 for Cs. In order to transport 
the Cs atoms by a defined distance, a precisely controlled relative detuning <5 between the lattice beams is used, which 
causes the standing wave interference pattern to move at a velocity v = Aiat^/2 for a specific amount of time |37| . For 
details also see Section inii 
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FIG. 2. Loading curve and histogram of the high gradient MOT. (a) Typical loading trace of the single atom MOT. 
After 8s a first atom is loaded into the MOT and in steps of 3-5s an atom is gained. The time is binned in 110 ms to collect 
enough fluorescence photons, (b) Histogram of the electron count rate taken over 30 traces with 267 data points each. The red 
solid line shows a Multi-Gaussian fit to the data, (c) Probabilities for each atom number obtained from the area of each peak 
in the multi-Gaussian fit (gray bars) and from a Poisson distribution (red bars) with mean atom number n=1.23, calculated 
from the histogram. 


Combining Single Atoms with the Quantum Gas 


During evaporation, the Rb cloud is prepared in the magnetic held insensitive state = 0) by a radio 

frequency transition, see Section m After it has been sufficiently cooled down to be well localized in the dipole trap 
crossing region, the Cs MOT is switched on. To avoid immediate losses of the single Cs atoms due to light-induced 
collisions [381440) . both species are trapped in different traps with a displacement of 110 ^m. Figure [ll| shows the 
probability of Cs atoms surviving a 1.3 s MOT period while a cold thermal Rb cloud is present in the dipole trap 
crossing region in dependence of the crossing region’s position. The minimum in survival probability marks the 
position of maximum overlap between MOT and dipole trap. 

As an initial step towards immersing the Cs atoms into the Rb BEC, we immerse them into a thermal Rb gas. Here, 
the horizontal dipole trap is sufficiently deep to directly trap Cs atoms from the MOT. After switching off the Cs 
MOT, we let the two species interact for a certain time before pushing Rb out of the trap by shining in resonant light 
perpendicular to both dipole trap beams. The high gradient MOT is switched on to recapture the Cs atoms from the 
dipole trap and the number of atoms in both MOT phases is compared to each other. Figure [IT] shows a measurement 
of the survival probability in dependence of the time between switching off the Cs MOT and pushing out the Rb, 
which we refer to as ’interaction time’. For a Rb cloud in |F = l,mF = 0) with 38 x 10^ atoms at a temperature of 
T = 3 /rK at trap frequencies of w = 27r x (63,1300,1300) Hz and a peak density of po = 1.3 x 10^^ cm“^, we measure 
a lifetime of r = 112 ms for Cs atoms in \F = 3). This is lower than expected because with a 3-body loss coefficient 
of L 3 = 5 X 10“^®cm®s“^ |24j the 3-body loss-rate can be calculated to be V'^s-body = L 3 x Pq = i/nsms at peak Rb 
density po and the Cs atoms should see a lower average density. We have not yet fully examined which processes limit 
the lifetime, but we expect two-body losses to play a role since the two species are not yet optically pumped to their 
lowest Zeeman substate. 
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Position of dipole trap crossing region [pm] 


FIG. 3. Light assisted collisions. The graph shows the probability of Cs atoms surviving a 1.3 s MOT period while a cold 
Rb cloud is present in the dipole trap crossing region in dependence of the crossing region’s position along the horizontal dipole 
trap beam. We attribute the losses to light assisted collisions of Cs with Rb in the presence of near resonant Cs MOT light. A 
negative Gaussian fit (red) serves as a guide to the eye. 



Interaction time [ms] 


FIG. 4. Lifetime of Cs in presence of a cold Rb cloud. The graph shows the probability of recapturing Cs atoms from 
the dipole trap after they have been stored together with a cold Rb cloud for a given ’interaction time’. Without Rb present in 
the trap, the survival probability is 0.69 (dashed gray). For a Rb cloud in [T = 1, mp = 0) with 38000 atoms at a temperature 
of T = 3 /iK at trap frequencies of oj = 27r x (63,1300,1300) Hz and a peak density of po = 1.3 x 10^^ cm“®, the fit reveals a 
lifetime of r = 112 ms for Cs atoms in [T = 3). 


III. METHODS 


Our experiments take place in a two-chamber vacuum system (see Figure III and III) consisting of a low pressure 
(« 10“^*^ mbar) and a high pressure region (« 10“^ mbar), separated by a differential pumping section of length 
83 mm and diameter increasing from 1.8 mm to 4 mm. The low pressure region is formed by a glass cell, whereas the 
high pressure region is located in a titanium chamber and contains the 2D MOT for ®^Rb which is loading atoms 
from the background gas. The distance between the 3D MOT region and the 2D MOT is approximately 30 cm. 


The coil system shown in Figure [Hi] b) is installed at the low pressure side, providing a quadrupole field for the 
Cs and Rb MOTs, a homogeneous field to address Feshbach resonances, and a vertical magnetic field gradient for 
Stern-Gerlach experiments. Additionally, compensation coils are installed to provide weak, homogeneous fields in all 
three dimensions up to 2 G to compensate magnetic stray fields. 
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FIG. 5. Schematical top view of our 2D/3D MOT system. 
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FIG. 6. Overview of the experimental setup, a) Experiment in top view: With dichroic mirrors beams for the optical 
dipole trap, MOTs for each species, the optical lattice for Cs, and an optical pumping beam for Rb are overlapped. A CCD 
camera resembles the position of the cold atoms which is useful for the initial alignment of the dipole trap and the optical 
lattice to the position of the Cs MOT. b) Front view with main magnetic coils (dimensions in mm), c) Side view. 


Cooling and Trapping 

Magneto Optical Traps for Rb The 2D MOT [41] setup is based on a design described in [42l|43|. The elongated 
titanium vacuum chamber features four windows with a clear aperture of 110 mm x 25 mm as an optical access for 
the MOT beams. Four distinct cooling regions are formed by orthogonal retro-reflected beams with diameters of 
2 X itiQ = 25 mm, and a viewport allows for an additional near-resonant beam (’push-beam’) along the cooling axis 
to increase the flux of the atomic beam. It is 1.9 MHz blue detuned from the F = 2 —>■ F' = 3 transition with an 
intensity of 120 The windows are surrounded by two quadrupole coil pairs, each coil consisting of 130 windings 

of 1 mm enameled copper wire. We typically operate the MOT at magnetic field gradients of 11 G/cm in vertical and 
horizontal direction. Our 3D MOT features a standard six beam configuration, where detuning and intensity can be 
controlled independently of the 2D MOT beams. We usually load the 3D MOT from the 2D MOT atomic beam for 
2.5s at a loading rate of 10® atoms/s before we enter a compressed MOT (’CMOT’) phase jUJHSj (see Figure [in| . The 
parameters for these Rb MOTs are listed in Table [^ 
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FIG. 7. CMOT scheme. We achieve an optimal transfer rate of Rb from the MOT to the dipole trap by increasing the 
red-detuning of the cooler while reducing the power of the repumper. 


cooler power 
cooler detuning 
repumper power 
repumper detuning 
beam size wo 
magnetic field gradient 


2D MOT 3D MOT CMOT 

240 mW 100 mW 50 mW 

-13.2 MHz -13.9 MHz -58 MHz 
28 mW 30 mW 14.8 pW 

-5 MHz -5 MHz -5 MHz 

12.5 mm 10 mm 10 mm 

11 G/ cm 6 cm radially 

11 G/ 

cm axially 


TABLE I. Parameters for the operation of different MOTs for ®^Rb. 


The laser light for operation of Rubidium MOTs, internal state manipulation and detection is provided by the 
master oscillator power amplifier (MOPA) setup shown in Figure [TIT} We employ interference-filter-stabilized extended- 
cavity diode lasers [46] (Radiant Dyes NarrowDiode), frequency stabilized to a frequency-modulation spectroscopy 
of Rubidium 87. The lasers’ output power of « 50 mW is amplified by home-made tapered amplifiers (TAs) up to 
2W, frequency shifted by acousto-optical modulators (AOMs) and transported to the experiment via polarization 
maintaining, optical single-mode fibers (Nufern PM780-HP). 

Single Atom MOT The position of the Cs MOT is overlapped with the Rb MOT as both use the same coil system. 
A set of 6 diaphragms with variable aperture mounted in the coil holders is used to align the MOT beams to precisely 
overlap in the middle of the glass cell. The cooling light is 10 MHz red detuned to the F=4 —)■ F’=5 transition of the Cs 
D2 line and has a total power of typically 500 /iW. The repumping light is on resonance to the F=3 —>■ F’=3 transition 
with a total power of typically 15/rW. In every beam pair a piezo driven mirror at 110 Hz frequency destroys phase 
coherence between orthogonal beam pairs and therefore avoids interference effects, which lead to an instable MOT 
position. To keep the MOT loading time short, a low magnetic field gradient of 40 G/cm axially and 20 G/cm radially is 
applied for 20 ms. During this sufficiently short time the trap volume is large enough to load on average one atom m- 
In a next step the magnetic field gradient is increased in 8 ms to 275G/cm in axial and 140G/cm in radial direction. 
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FIG. 8. Schematic diagram of our laser setup for Rubidium trapping, manipulation and detection. The notation 
F ^ F' refers to the hyperfine transition of the D 2 line to which the light is resonant, while F —>■ {F[, F 2 ) refers to a crossover 
signal of the Doppler-free spectroscopy. 



FIG. 9. Evaporation scheme. The graph shows the horizontal beam’s power (black) and the vertical beam’s power (orange) 
over time, including parameters of the exponential ramps. 


This effectively pins the atom number and avoids additional atom loading during the imaging process. The duration 
of the low gradient stage depends on the vapor pressure of the Cs atoms and the required atom number. For an 
efficient transfer of the trapped atoms into the optical dipole trap, a low MOT temperature is needed. Therefore we 
increase the red detuning of the cooler light to 72 MHz for 50 ms, while setting its power so low that we just do not 
lose the atoms. We release the atoms into the dipole trap by switching off the MOT beams. The cooler is switched 
off 2 ms after the repumper ensuring that the atoms remain in their lowest fine-structure state \F = 3). 

Optical Dipole Trap Our optical dipole trap is a crossed beam trap at 1064 nm. The trap is formed by a horizontal 
beam with a focal waist of wq = 22 /rm at 4 W and a vertical beam with a focal waist of wq = 165 /rm at 12 W of power. 
We measure trap frequencies of 3.7 kHz radial to the horizontal beam and 120 Hz along the horizontal beam, while 
the axial confinement is mainly caused by the vertical beam. The beam setup allows both for a forced evaporation 
scheme [15] as well as standard, passive evaporation [15]. We keep the pointing instability minimal by using optical 
fibers to transport the light to the experiment (LMA-PM-15 by NKT Photonics for the vertical beam and Liekki 
Passive-lO/125-PM for the horizontal beam). AOMs between the laser light source (Nufern NUA-1064-PD-0050-D0 
fiber amplifier) and the optical fibers are controlled by a PID controller that determines and stabilizes the laser power 
at the experiment with a bandwidth of 110 kHz. 

Evaporative Cooling Both beams are switched on at full power during the CMOT phase. At the end of the CMOT 
phase, the MOT lasers and the quadrupole field are switched off. The repumping light is switched off 200 /rs before 
the cooling light, so that the atoms are pumped to the F = 1 state and approx. 10® atoms are transferred to the 
dipole trap. Figure |ITT| shows our subsequent evaporation scheme. After 300 ms of self-evaporation we decrease the 
horizontal beam’s power exponentially in two sections with different decay constants r. As soon as the increase in 
density within the crossing region leads to three-body losses, we also ramp down the power of the vertical beam 
until a BEG with 2.5 x 10^ atoms forms after a total evaporation time of 4.3 s. The entire BEG sequence, including 
3D MOT loading, GMOT phase and evaporation ramps has been optimized with an evolutionary algorithm, that is 
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FIG. 10. Driving a Landau-Zener transition and performing a Stern-Gerlach experiment. Rb atoms in the 
\F = l,mF = 1) state are transfered to the \F — l,m_F = 0) state by performing a Landau-Zener sweep. The magnetic field is 
ramped such that the detuning of the radio frequency with respect to the Imp = 1) —>■ |mF =0) transition changes from far 
blue detuning to far red detuning. After the atoms are released from the dipole trap, a magnetic field gradient is applied in 
order to perform a Stern-Gerlach Experiment. 


implemented as a part of our timing software m- 

Internal State Preparation The Rb atoms are prepared in the magnetically insensitive state 1^=1, mp = 0) during 
evaporation so that the magnetic field gradient of the single atom MOT does neither heat nor destroy the cold Rb 
cloud: During self-evaporation, a magnetic field of 1.4 G is switched on to lift the degeneracy of the Zeeman substates. 
Then the atoms are pumped to the |F = l^mp = 1) state by a 900/rs light pulse resonant to the F = 1 —)■ F' = 1 
cr■'■-transition and a 500 ps light pulse resonant to the F = 2 —)■ F' = 2 7r-transition of the D 2 line. Both beams are 
turned on simultaneously. Their respective powers are 300 nW and 480/xW with red detunings of 2.6 MHz and 1 MHz 
at equal beam waists of 1.2 mm. 

By applying a magnetic field gradient of typically during a 13 ms time of flight experiment (see Fig. HI I, 

we perform a Stern-Gerlach experiment to measure the population of magnetic substates. Without optical pumping 
we observe an almost equally distributed spin mixture, while with optical pumping approximately 90 % of the atoms 
are In the mp = 1 Zeeman substate and approximately 10% remain in mp = 0. Optical pumping does not deplete 
the number of atoms in the condensate. 

After optical pumping we transfer the population to the mp = 0 state with a Landau-Zener sweep as shown in 
Figure [ITT] We apply a radio frequency at 10.176 MHz for 100 ms while we increase the magnetic field linearly. The 
magnetic field ramp is chosen such that the detuning of the radio frequency with respect to the \mp = !)—:► \mp = 0) 
transition falls from approximately -1-550 kHz to —550 kHz. We observe a Rabi-frequency of H = 27r x 2.1kHz with a 
transfer efficiency of 95%, with no negative impact on the number of atoms in the BEG. 


Imaging 


We image the Rb cloud by absorption Imaging with linearly polarized light resonant to the |F = 1) —)• |F' = 1) 
transition of the D 2 line, see Figures [HT] and |lTT] 

To observe the single atoms we use fluorescence Imaging. Near resonant light from the MOT beams is scattered 
in the whole solid angle. We collect 3.3% of the photons corresponding to 0.012 pW/atom for a saturation parameter 
= ///s = 1 by a custom made objective. It has a numerical aperture (NA) of 0.36 and is placed beneath the glass 
cell at a distance of 30.3mm to the atoms’ position (see Figure HI:). Stray light protection is crucial and the whole 


beam path Is located Inside black anodized lens tubes and mirror housings. 

The objective tube is mounted on a precision xyz linear stage (Newport Ultralign Model 562-XYZ) for accurate 
positioning. The light which is collimated by the objective is focused onto an EMCCD camera (Andor iXon 3 897) by 
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.—fl—^—fl—. .- f, -^-f, 

fi = 120 mm ^2 — 300 mm 


FIG. 11. Absorption imaging with a divergent beam. The drawing shows both virtual rays originating from the 
atomic cloud and real rays from the divergent imaging beam. The atomic cloud lies in the focal plane of the imaging system. 
Independent of where the divergent beam’s focus is located, a shadow magnified by cast onto the CCD chip. 



FIG. 12. Fluorescence images of single Cs atoms in MOT and lattice. Four consecutive images of three single Cs 
atoms, a) Fluorescence image of Cs atoms trapped in a MOT, b) after being transferred to the species-selective conveyor belt 
lattice, c) after being moved for 30 pm, d) after being held in place, e) after being transferred back to the MOT. 


a lens with focal length / = 1000 mm. This yields a magnification of M=33. The camera chip has 512 x 512 pixel with 
a pixel size of 16 x 16 pm which leads to a field of view of 250 pm. An advantage of this setup is that the objective 
can be aligned independently from the focusing lens which is producing the image on the camera. The high EMCCD 
gain and the quantum efficien cy o f nearly 60% of the camera allows to observe single atoms in the MOT as well as in 
the optical lattice (see Figure [TTT| ). During imaging the atoms are illuminated by the repumper light of the MOT and 
cooler light, which is 20 MHz red detuned in order to compensate the light shift caused by the dipole trap. Exposure 
times of 150 ms (400 ms) in the MOT (lattice) are usually used to image the atoms. The number of atoms in the 
MOT can be determined with close to 100 % reliability from the brightness of the image. 


Optical Lattice 

Our optical lattice consists of a two beam setup (see Figure [nT| ), where both beams are guided to the experiment by 
means of optical fibres with identical length. This helps to reduce phase drifts and guarantees a good beam quality at 
the experimental side. To make sure that both beams hit the small-volume Cs MOT, we exploit the AC stark effect. 
The dipole trap is red detuned and hence reduces the fl uore scence light of the MOT, whereas the lattice frequency is 
blue detuned and enhances the fluorescence (see Figure [hT]). 


Lattice Transport 

For the single atom transport in the conveyor belt optical lattice m, a relative detuning S between the two lattice 
beams is employed to create a standing wave pattern, moving at velocity v = Aiat<J/2. For a typical transport, the 
detuning between the beams is ramped in a trapezoid shape: a linear ramp from zero to the maximum detuning 
yields an acceleration of the atoms in the lattice, followed by a plateau of constant detuning where the atoms move 
at constant velocity. After this plateau, the detuning is ramped back to zero and the atoms are decelerated again. 
The transport distance is given by the integral over the velocity. Thereby large transport distances in the range 
of millimeters can be realized, which are only limited by trap size. This is in contrast to phase shifting transport 
approaches, where the maximum transport distance is given by the maximum phase shift EH. The absolute maximum 
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FIG. 13. AC Stark effect of dipole trap and lattice. A typical fluorescence trace with lattice and dipole trap beam. 
During the MOT phase one atom is captured. Switching on the dipole trap reduces the fluorescence level by about 70%. In 
combination with the lattice a slight increase is observed because it compensates the light shift of the dipole trap. For the 
combination of MOT and lattice an increase of about 10% in fluorescence level is observed. In the end the count rate due to 
stray light and the dark count rate is detected. A binning time of 100 ms is used. 



FIG. 14. Lattice control setup. The DDS-based signal generation (blue boxes) supplies two phase stable radio-frequencies 
to the amplifier chains (red boxes) which drive the lattice AOMs. 


acceleration Omax = kUo/mcs ~ 4 • 10® m/s^ in the lattice is determined by the Cs mass tocsj the laser wave number 
k and the potential depth Uq « 1.5 mK of the standing wave potential. 

In order to control the detuning, each lattice beam is frequency shifted by a common value of /i ~ /2 = 160 MHz 
in an AOM double-pass setup. The relative detuning between the two beams then is given by <5 = /i — /2 and can 
be controlled by the two RF frequencies supplied to the AOMs. Both RF signals have to be phase-stable compared 
to each other and the offset-frequency of both channels has to be equal up to a fraction of a Hertz to yield a stable 
standing wave pattern when both beams interfere for 6 = OHz. Therefore, we use a driver electronics based on direct 
digital synthesis (’DBS’) with an amplifier chain shown in Figure HI The DBS chips employed (AB9954) sample 
the output sine wave in a digital circuit and have an analog digital converter stage for signal output. By supplying 
both BBS chips with the same clock signal which is locked to a Rb frequency standard, the relative phase stability 
between the signals is guaranteed by the digital sampling of the output signal. The frequency of the output sine wave 
is set by a 32 bit control parameter, yielding a frequency resolution of 0.09 Hz which allows for very small frequency 
variations compared to the output frequency of 80 MHz. To supply output powers of up to ^30 dBm a two stage 
amplifier chain combined with a voltage controlled attenuator for output power control is applied. The RF power 
level can be directly used to adjust and stabilize the light intensity of the lattice beams. The BBS chips are controlled 
by a microcontroller which provides exact timing of frequency changes and stores the data for the frequency ramps. 


IV. CONCLUSION 

Our experimental results constitute an ideal starting point for studying coherent interactions of single Cs impurities 
in a Rb BEC. We have shown that, for this system, the experimental realization, control and detection of a Bose 
polaron in the weak and intermediate coupling regime of the Frohlich model is accessible. Moreover, our experimental 
apparatus not only facilitates the creation of a tunable Bose polaron but it also allows for controlled positioning and 
dynamical studies of impurities in a quantum fluid. 
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